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ABSTRACT 

We calculate the colours and luminosities of redshift z = 0.1 galaxies from the EA¬ 
GLE simulation suite using the GALAXEV population synthesis models. We take into 
account obscuration by dust in birth clouds and diffuse ISM using a two-component 
screen model, following the prescription of Chariot and Fall. We compare models in 
which the dust optical depth is constant to models where it depends on gas metallicity, 
gas fraction and orientation. The colours of eagle galaxies for the more sophisticated 
models are in broad agreement with those of observed galaxies. In particular, eagle 
produces a red sequence of passive galaxies and a blue cloud of star forming galaxies, 
with approximately the correct fraction of galaxies in each population and with g — r 
colours within 0.1 magnitudes of those observed. Luminosity functions from UV to 
NIR wavelengths differ from observations at a level comparable to systematic shifts 
resulting from a choice between Petrosian and Kron photometric apertures. Despite 
the generally good agreement there are clear discrepancies with observations. The blue 
cloud of EAGLE galaxies extends to somewhat higher luminosities than in the data, 
consistent with the modest underestimate of the passive fraction in massive eagle 
galaxies. There is also a moderate excess of bright blue galaxies compared to obser¬ 
vations. The overall level of agreement with the observed colour distribution suggests 
that EAGLE galaxies at z = 0.1 have ages, metallicities and levels of obscuration that 
are comparable to those of observed galaxies. 

Key words: galaxies: scaling relations, galaxies: simulation, galaxies: luminosity 
function, galaxies: colour-magnitude diagram 


1 INTRODUCTION 

The basic scenario for how galaxies form and evolve is well 
established: gas accretes onto deepening dark matter po¬ 
tential wells, cools and makes stars. Although this basic 


paradigm has been accepted for many years (e.g. Rees & Os- 


triker|1977 White & Rees 1978 White fc Frenk|1991[ ), many 
important aspects are still poorly understood. For example, 
the shape of the galaxy stellar mass function (GSMF) and 
of the halo dark matter mass function are quite different, so 
a simple formation model in which a halo of given mass con¬ 


tains a galaxy whose stellar mass is equal to a fixed fraction 
of its halo mass is ruled out (e.g. [White fc Frenk|199l| |Ben-| 
son et al. 20031. In addition, galaxy surveys have revealed 


the presence of many correlations between the stellar prop¬ 
erties of galaxies and galaxy mass, for example bimodality 
in the colour-magnitude diagram in the form of a blue cloud 
of star forming galaxies at lower mass, and a red sequence 
of mostly passive galaxies at higher mass (e.g. [Baldry et al.j 
20041. Clearly, the dynamics and interaction physics of gas 


play important roles in determining the properties of galax- 
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It is thought that at the faint end of the GSMF super 


novae (SNe) quench star formation in small galaxies (Larson 
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1974 Dekel & Silk 19861 with reionisation effectively pre¬ 


venting galaxies from forming in dark matter halos below a 


minimum mass (e.g. Rees 1986 

Efstathiou 1992 

Thoul & 

Weinberg|[1995 Okamoto et al. 

2008). Together these pro- 


cesses shape the GSMF at low stellar masses, with dark 
matter halos containing increasingly feeble galaxies with de¬ 


creasing halo mass, until most halos remain dark (Sawala 


et al. 20141. In contrast, at higher galaxy masses it is thought 


to be the feedback from accreting black holes (BHs) that in¬ 
troduces a near exponential cut-off in the GSMF (|Bowei 
et al.|2006 Groton et al.|2006 |. 


That feedback from star formation and BHs together 
shape the galaxy stellar mass function is energetically plau¬ 
sible, and semi-analytical models (e.g. Cole et al.|2000 Hen- 


riques et al.||201^ |Gonzalez-Perez et al.||2014| [Porter et al. 

2014) and numerical simulations (e.g. |Oppenheimer et al. 
2010 [Puchwein fc Springel||2013| [Vogelsberger et al.||2014[ 
Schaye et al.poTs ) that appeal to these processes are able 
to produce model GSMFs that compare well to observations 
at redshift 2 ~ 0. However, neither the semi-analytical mod¬ 
els nor these simulations come anywhere near resolving the 
scales at which SNe and BHs inject energy, and so cannot 
a priori compute the net efhciency of the resulting feed¬ 
back. Simulations therefore need to rely on a phenomeno¬ 
logical description of these crucial processes occurring on 
unresolved (‘subgrid’) scales, using observations to calibrate 
the parameters that appear in the subgrid modules. It is 
then important to quantify the uniqueness and degeneracies 


in such modelling (Schaye et al. 2010 Crain et al. 2015), 
while at the same time use very high-resolution simulations 


(e.g. Hopkins et al.|2011[ Creasey et al.|2013 2015 Martizzi 


et al. 2014 Rosdahl et al. 2015) to try to bridge the gap 


between numerically unresolved and resolved scales. 

The Evolution and Assembly of GaLaxies and their 
Environment (eagle [Schaye et al.|[20T^ [Crain et al.||201^ 
suite of simulations uses the 2 ~ 0.1 GSMF, together with 
observations of galaxy sizes, to calibrate the subgrid physics 
modules. The GSMF is however not directly observable, but 
is inferred from luminosity functions after applying correc¬ 
tions for dust obscuration, and using simple stellar popula¬ 
tion (SSP) synthesis models that involve assumptions about 
stellar evolution, star formation histories, metallicity depen¬ 
dence of stellar emission, etc. With simulations such as EA¬ 
GLE we can take the converse approach, attempting to re¬ 
produce the observational relations by inputting physical 
quantities tracked by the simulation. This has the advantage 
of allowing one to use properties modelled self-consistently 
such as gas content, metallicity and age to derive observable 
quantities, rather than treating them as free parameters in 
‘SED fitting’ (e.g. Walcher et al.|2011 1 to estimate physical 
properties from observations. The colours of EAGLE galax¬ 
ies are also an important test of the realism of the fiducial 
EAGLE model. 

In this paper we examine to what extent mock lumi¬ 
nosities computed from EAGLE galaxies using SSP models 
and a simple correction for dust obscuration reproduce the 
observed luminosity functions (in a range of broad bands), 
as well as galaxy colours. The aim is twofold: to provide a 
test of the realism of EAGLE, but also to test at some level 
whether the procedure of going from luminosity to stellar 
mass is reliable, as investigated in a recent paper by [Torre3^ 
et al. (20141 using SED fitting of galaxies from the Illustris 


Table 1. Numerical parameters of those simulations of the EA¬ 
GLE suite that are used in this paper. Prom left to right: sim¬ 
ulation identifier, side length of cubic volume L in co-moving 
Mpc (cMpc), initial mass mg of baryonic particles, Plummer- 
equivalent gravitational softening tprop at redshift 2 = 0 in proper 
kpc (pkpc). 


Name 

L 

mg 

£prop(^ — 0) 


cMpc 

M© 

pkpc 

Ref-L025N0376 (Ref-25) 

25 

1.81 X 10® 

0.70 

Recal-L025N0752 (Recal-25) 

25 

2.26 X 10® 

0.35 

Ref-L100N1504 (Ref-100) 

100 

1.81 X 10® 

0.70 


simulation | Vogelsberger et al.|2014 1. It is perfectly possible 
that EAGLE galaxies have the wrong colours but the right 
stellar masses and stellar ages if errors in dust modelling 
are severe. However, if masses and colours agree with the 
data, then this increases our confidence that we can use EA¬ 
GLE to investigate, for example, the origin of the observed 
bimodality of galaxy colours, or the nature of the galaxies 
that lie in between the red sequence and blue cloud in a 
colour-magnitude diagram. 

In what follows we compare to photometric data from 


the Galaxy And Mass Assembly survey (gama Driver et al. 


20111. This is a spectroscopic follow-up based on photomet¬ 


ric data from the Sloan Digital Sky Survey (SDSS) and the 
UK Infrared Digital Sky Survey (UKIDSS), with details on 
the targeting and star-galaxy separation in [Baldry et al.j 
(2010) and on the GAMA-processed photometry, including 
matched aperture photometry from u to AT in [Hill et ah] 
(20111. The GAMA survey has been designed for high uniform 
spectroscopic completeness ( Robotham et al.|2010 1 and pro¬ 
vides accurate redshifts (using AUTOz; [Baldry et al.[[2014 ) 
for a catalogue of ~ 190, 000 galaxies, as presented byjTaylorj 


et al. (2014). 


This paper is organised as follows. We begin with an 
overview of the EAGLE simulations, with emphasis on those 
aspects that are most relevant for the SSP modelling. In !|^ 
we detail the development of our photometric model, con¬ 
centrating on emission and absorption in §3.1| & §3.2| re¬ 
spectively. This model is applied to yield an optical colour- 
magnitude diagram (CMD) and multi-band luminosity func¬ 
tions (LFs) for galaxies, which are plotted and discussed in 
[Q We discuss our findings in ^and summarise in !|^ 


2 THE EAGLE SIMULATIONS 


Full details of the eagle simulations can be found in jSchayej 
[et al.[ | [2015[ | and [Crain et al.[ ( [2015[ ) (hereafter S15 and 
C15 respectively); here we give only a brief overview. The 
EAGLE simulation comprises a suite of cosmological hydro- 
dynamical simulations of periodic cubic volumes performed 
with the Gadget-3 TreeSPH code (which is an update of 
the Gadget-2 code last described by Springel||2005 1. Sim¬ 
ulations were performed for a range of volumes and nu¬ 
merical resolutions. Here we concentrate on the reference 
model, using simulations at different resolution to quan¬ 
tify numerical convergence. The reference model assumes 
a ACDM cosmology with parameters derived from the ini- 
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tial Planck ( Planck Collaboration et al.|2014 l satellite data 
release (fib = 0.0482, fidark = 0.2588, JIa = 0.693 and 
h = 0.6777, where Ho = 100 h km s“^ Mpc“^). Relevant 
properties are listed in Table 

We modified the treatment of hydrodynamics in 
Gadget-3 to use the conservative pressure-entropy SPH 


formulation of Hopkins (20131, the artificial viscosity switch 
introduced by Cullen fc Dehnen (|2010(, an artificial con¬ 


duction switch inspired by Price (20081, the kernel from 


Wendland (19951 and the timestep limiter of Durier & Dalla 


Vecchia (20121. Motivation and tests of this ‘Anarchy’ ver¬ 


sion of SPH are presented in Dalla Vecchia (in preparation, 
see also Appendix A of S15), whereas Schaller et al (in prepa¬ 
ration) examine the effects of using this modified version of 
SPH on galaxy properties. 

A crucial aspect of EAGLE is that the parameters de¬ 
scribing the subgrid modules have been calibrated on the 
observed 2 ~ 0 GSMF and galaxy sizes. This good agree¬ 
ment extends to many other observables that were not con¬ 
sidered during the calibration, such as specific star formation 
rates (S15), the evolution of the GSMF ( Furlong et al.|20l4 1, 
molecular hydrogen fractions (Lagos et al. 20151, and ab¬ 


sorption by intergalactic metals and neutral hydrogen (S15 
Rahmati et al.|[2015 |. 

We now briefly describe those aspects of the subgrid 
model most relevant for this paper. 


motivated in Dalla Vecchia & Schaye (20121, and adapted 


for the simulation as described in S15 and C15. The forma¬ 
tion of supermassive black holes and gas accretion onto them 
is implemented as in Booth & Schaye (20091, with modifi¬ 


cations described in Rosas-Guevara et al. (20131 and S15, 


using a single feedback mode. 


2.2 Identifying galaxies 


To group star particles into ‘galaxies’, we proceed as fol¬ 
lows. We begin by identifying dark matter halos using the 
friends-of-friends algorithm (Davis et al.||198^( with a link¬ 


ing length of 0.2 times the mean dark matter inter-particle 
spacing to identify regions that are overdense by a factor 


of ~ 200 (Lacey fc Cole 

1994|). We then use the subfind 

algorithm (|Springel et al. 

|2001| |Dolag et al.||2009| to iden- 


tify self-bound substructures (subhalos) within halos of dark 
matter, stars and gas, which we identify with galaxies. Mas¬ 
sive galaxies in EAGLE have extended density profiles. To 
assign luminosities to mock galaxies, we only include light 
emitted within a sphere of radius 30 proper kpc (pkpc), 
centred on the minimum of potential of each subhalo. The 
motivation for the choice of aperture is discussed further in 
section TS.l.dl and also in S15. 


2.1 EAGLE Subgrid physics 

The subgrid modules of EAGLE are partly inspired by the 
OWLS and GIMIC simulations ( |Schaye et al.|20ro||Crain et al. 
2009). Star formation is implemented as described in Schaye 


fc Dalla Vecchia (2008): above the metallicity-dependent 
star formation threshold of |Schaye| ( |2004[ ), cold enough parti¬ 
cles (T ~ 10‘^K) are converted to star particles stochastically 
at a pressure-dependent rate that reproduces the observed 
Kennicutt-Schmidt star formation law. Each simulation star 
particle is assumed to represent a coeval population of stars 
formed with a Chabrier] ( |2003[ ) IMF, comprising stars with 
masses in the range [0.1,100] Mq. 

Stellar evolution is implemented as described in S15 and 
Wiersma et al. ( 2009a[ ). We follow the production and re¬ 
lease into the interstellar medium of 11 elements from three 
channels of stellar evolution (AGB stars, type I and type H 
SNe) using metallicity-dependent stellar lifetimes and stel¬ 
lar yields. We also track a separate ‘total’ metallicity (the 
mass fraction of elements more massive than helium), Z, 
to account for elements not tracked explicitly. When a gas 
particle is converted into a star particle it inherits the gas 
particle’s abundances. In addition to a particle metallicity, 
the simulation tracks smoothed metallicities which are com¬ 
puted using the SPH formalism to partly remedy the ab¬ 
sence of mixing in the calculation (see Wiersma et al.|2009a 
for motivation). The ages, masses and metallicities of the 
star particles are the main ingredients of the simple stellar 
population (SSPs) models used below. 

Radiative cooling and photo-heating in the presence of 
an optically thin UV/X-ray background, as computed by 
Haardt fc Madau (|2001| ), is accounted for as described in 
Wiersma et al.| ( |2009b[ ). The crucial processes of feedback 
from star formation is implemented by stochastically heating 
particles by a fixed temperature increment, as described and 


3 PHOTOMETRY 

This section explains how we compute luminosities and 
colours for the simulated EAGLE galaxies. We begin by mod¬ 
elling luminosities of star particles, taking into account their 
ages and metallicities using population synthesis (section 
3.1), the photometric system used to calculate magnitudes 


for direct comparison to observation (section 3.1.3), and 


the effects of dust absorption (section |3.2[ ). The results of 
this section are summarised in Fig. in which we plot his¬ 
tograms of g — r colours of EAGLE galaxies in narrow stellar 
mass bins for different models, ranging from a simple emis¬ 
sion model without dust, to a model including a multivariate 
treatment for dust. 


3.1 Source Modelling 

Below we compute luminosities of EAGLE galaxies from the 
ultraviolet (UV) to the near infrared (NIR). We limit our¬ 
selves to modelling stars, neglecting both nebular emission 
lines and light from AGN. Light absorbed by dust is as¬ 
sumed to be re-emitted in the far infrared which we do not 
study in this paper. As we also neglect scattering by dust, 
we treat individual wavelength bands independently. This 
approximation is supported by observations showing that 
scattering is a negligible contributor to the observed atten¬ 
uation curve in galaxies (eg. Fischera et al.|2003 ). 

Population synthesis models provide spectra for a dis¬ 


crete range of simple stellar populations (SSPs) (e.g. Bruzual 
fc Chariot 2003 Maraston 2005 Leitherer et al.|1999 ). These 
SSPs represent populations of stars characterised by their 
total initial mass, formation time, and composition while 
assuming some stellar IMF. By decomposing the stellar com¬ 
ponent of a galaxy into a superposition of SSPs, the spectral 
energy distribution of an entire galaxy can be approximated. 
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4 Trayford et al. 


We treat each eagle star particle as an SSP with given ini¬ 
tial stellar mass, metallicity and age. 


3.1.1 SSP Ingredients 


Given our implementation of star formation, where gas par¬ 
ticles are wholly converted into star particles, the typical 
mass of an EAGLE star particle is ~ lO^M©. Stars are as¬ 
sumed to form with a Chabrier ( 2003| IMF (for consistency 
with the evolutionary models used in eagle), and they in¬ 
herit the SPH-smoothed metallicity, Z, from their parent 
gas particle. The mass of a newly formed star particle is 
purely set by numerical resolution; the particle should not 
be thought of as representing a star cluster. In fact, 10® 
Mq is higher than the stellar mass formed in giant Hll re¬ 
gions (e.g.|Relano fc Kennicutt|2009[|Zaragoza-Cardiel et al.| 


2014). This poor sampling of star formation can adversely af¬ 


fect luminosities of EAGLE galaxies. Indeed, a single recently 
formed star particle can affect the colour of a galaxy. We 
try to mitigate this numerical artifact by employing a finer 
sampling of the star formation history of recently formed 
stars, as described in Appendix We note that this has 
only a small effect for optical colours and thus for the re¬ 
sults presented here. 

The metallicity of stars affects their colours resulting 


in the well-known age-metallicity degeneracy (e.g. Worthey 


1994). In addition, Z affects stellar evolution leading to dif¬ 


ferences between models, particularly for the AGB phase 


(e.g Inoue 2012 Stancliffe & Jeffery 2007). In addition, 
metallicity of stars in EAGLE galaxies is set by the intricate 
interaction between enrichment of the ISM, gas accretion, 
and the extent to which galactic winds transport metals into 
the galaxy’s circum- and intergalactic medium. The details 
of such metal mixing are still poorly understood and nu¬ 
merically challenging to model. Hence, there is no a priori 
guarantee that EAGLE yields realistic stellar metallicities. 

The stellar mass - metallicity (M* — Z*) relations pro¬ 
vide a useful way of characterising stellar abundances as a 
function of galaxy mass, and as shown in S15, the Ref-100 
model yields stellar and gas-phase metallicities consistent 
with observations (Tremonti et al.|2004 Gallazzi et al.|20'05 


Zahid et al. 2014) for stellar masses M* > lO"'^^©. However, 


lower-mass galaxies in EAGLE tend to be more metal-rich 
than observed, with numerical resolution playing a role: the 
higher-resolution Recal-25 simulation agrees with the data 
for M* < lO^M©. It should also be noted that there are large 
systematics on the observed mass-metallicity relations (e.g. 
|Kewley fc Ellison||2008 ). We investigate the impact of stel¬ 
lar metallicity (Z*) on EAGLE colours in more detail in Ap- 
pendixj^ In our analysis, we use the EAGLE SPH-smoothed 
metallicities ( Wiersma et al.|2009b l for each particle, which 
yield less noisy estimates of Z^,. 


3.1.2 Stellar Population Synthesis (SPS) Modelling 

population synthesis models of 
, which provide the spectral en¬ 
ergy distribution (SED) per unit initial stellar mass of a 
SSP for a discrete grid of ages (ranging from t = 10® to 
2 X 10^® yr) and metallicities (ranging from Z* = 10“^ to 
0.05). We compute the SED for each stellar particle by in¬ 
terpolating the GALAXEV tracks logarithmically in age and 


We adopt the galaxev 
Bruzual fc Chariot] (2003 


Z*, and multiplying by the initial stellar mass. Values of 
Z* in the simulation span a wider range than the models of 


Bruzual fc Chariot (2003) encompass, and we extrapolate 
the model to predict colours for such extreme metallicitie^ 
We note that the more conservative approach of not extrap¬ 
olating introduces a bias. 

The GALAXEV spectra are based on stellar emission 
alone. These models are widely used, and have been shown 
to fit the local galaxy population in the optical bands with 
reasonable star formation and enrichment histories when 


used in conjunction with a dust model (e.g. Chariot & Fall 
2000 Cole et ar]|2000l. The choice of population synthesis 


model has been shown to be largely unimportant for low- 
redshift galaxy populations, especially in optical bands (e.g. 


Gonzalez-Perez et al. 2014). The effect of different models 
(e.g. |Maraston||2005[ [Conroy et ^|2009p , and in particular 
the more uncertain impact of thermally-pulsing AGB stars 
(TP-AGB) on the SED, should however be considered when 
surveying higher redshift (z > 1) galaxy populations (e.g. 
Maraston|2005 Gonzalez-Per^ et al.|20l4). 


The Bruzual & Chariot ( |2003| ) models specify Z* val¬ 
ues as absolute metal-mass fractions, where Z* affects the 
colours of stars through its impact on stellar structure and 
evolution - for example via the opacity, equation of state and 
mean molecular weight - and on stellar atmosphere models. 
Even so, the metallicity of the Sun (Z©) does enter the pop¬ 
ulation synthesis models because some parameters, such as 
the mixing length, are calibrated based on solar observables 
(Bressan et al. 1993). Recent literature determinations of 
Z© have shown significant variability, with a minimum of 
Z© ~ 0.0126 (e.g. Asplund et al.|2004 ) from the traditional 
value of Z© = 0.02 assumed in Bruzual & Chariot (2003). 
Although the eagle simulations do not make use of any 
particular solar abundance pattern or Z© value, a relatively 
low value of Z© = 0.0127 (Allende Prieto et al. 2001) has 
been assumed in analysis for consistency with |Wiersma et ah] 
( 2009b| ) (S15). The variation of Z© generally results from a 
different determination of the abundance of some important 
element, such as O, N, C or Fe which also implies a varia¬ 
tion of the abundance partition in the Solar model. In prin¬ 
ciple one should use galaxev SSP models with an abun¬ 
dance partition consistent with the assumed value of Z©, 
and take into account effects arising from the different mix¬ 
ing length calibration, to compute colours self-consistently. 
For now we neglect any such changes and use the original 
GALAXEV SEDs, as the effects of this change on broad-band 
colours are expected to be small, as long as one makes use of 
the absolute value of the metallicity (Bressan 2014 private 
comm.j^ 


3.1.3 Photometric System 

Given the SED for each star particle in the simulation, and 
a model for attenuation due to dust as a function of wave¬ 
length, we could compute the SED for each galaxy, and cal- 


^ We find that the effect of extrapolating metallicities, as opposed 
to clipping metallicities to that of the highest metallicity GALAXEV 
spectra, has a negligible effect on our results. 

^ We are grateful to S. Chariot and A. Bressan for their detailed 
explanation of the impact of Z© on the GALAXEV model. 
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culate a broad-band magnitude by convolving with a broad¬ 
band filter. Here we use a much simpler method, namely 
we first convolve the GALAXEV spectra (for each value of 
age and metallicity) with broad-band filters to obtain ‘un¬ 
obscured’ broad-band luminosities. We use these to obtain 
a broad-band luminosity for an EAGLE galaxy. We only then 
take into account dust attenuation (as described below). If, 
as we assume, the wavelength dependence of the dust at¬ 
tenuation is not very strong (he. the optical depth does not 
vary strongly over the wavelength extent of the filter), then 
these two approaches yield very similar results. For the dust 
models discussed below, we verified that this is indeed the 
case (section |3.2[ ). 

We use the ugrizYJHK photometric system for optical 
and near infrared photometry, to enable a direct compari¬ 
son to the GAMA survey (described in Driver et al.|2011 Hill 
et al.|20rT l. This survey is based on the photometry of SDSS 
(technical description in York et al. 20001 and UKIDSS, 


(technical description in Lawrence 20071. Note that when 


calculating photometry below, filter transmission curves in¬ 
clude atmospheric absorption to enable a like-for-like com¬ 
parison of simulation and observation. All magnitudes are 


calculated as rest-frame and absolute in the AB-system (Oke 


1974) in which the apparent magnitude mAB is defined as 


rriAB = -2.5 logio(F„) - 48.6 , 


( 1 ) 


where is the isophotal flux density (in cgs units) in a 
particular band (e.g. Tokunaga &: Vacca|2005|. 


3 . 1.4 Choice of aperture 

Individual EAGLE ‘galaxies’ are identihed as described in |2.2[ 
We select galaxies with at least 100 star particles, whose stel¬ 
lar mass is reasonably well converged numerically (S15). The 
line of sight is chosen consistently to lie along an axis of our 
simulation coordinates, yielding an essentially randomised 
orientation for each galaxy. 

Massive galaxies (M* > 10^°'®M©) in EAGLE have up 
to ~ 40% of their stellar mass in an extended halo beyond 
30 pkpc of the galaxy centre (‘intra-cluster light’ since most 
of these massive galaxies are at the centre of a group or 
cluster). Observationally, such galaxies also tend to have 
extended light distributions and, unsurprisingly, the lumi¬ 
nosity assigned to them depends on how such light is taken 


into account (e.g. Bernardi et al. 2013). We apply a con¬ 


stant aperture of 30 pkpc centred on the minimum of the 
gravitational potential of each subhalo for measuring the to¬ 
tal luminosity, L, of a galaxy. The luminosity and colour of 
a galaxy with a significant intra-cluster light component do 
depend on whether we apply a 2D aperture or a 3D aperture. 
This change in colour is due to colour gradients but also due 
to the inclusion or exclusion of small blue satellite objects 
below the signihcance of those identified by subfind. This is 
a similar issue to that encountered when isolating galaxies 
in astronomical data using software such as Sextractor 


(Jackson et al. 2010). We choose to apply a 3D spherical 
aperture, consistent with our previous analysis (S15, |Fur- 
|long et al.| ( |2014[ )). Such an aperture is shown in S15 to 
yield similar stellar masses to a Petrosian aperture, often 
used in observational studies. 

The aperture definition is not standardised in observa¬ 
tions, and can make a difference when a considerable frac¬ 


tion of stellar material lies outside the aperture. This is il¬ 
lustrated for the Kron and Petrosian apertures in |Driver| 
(2012), where luminous galaxies with high Sersic indices 


yield different magnitudes. Similarly, when analysing our 
simulations, the luminosities of EAGLE galaxies with M^, > 
lO^^M© are sensitive to the exact choice of aperture size. 
However, this is not the case for lower mass galaxies, for 
which the fraction of light in an extended halo is much lower. 


3.1.5 Model N 

The procedure for obtaining EAGLE galaxy photometry out¬ 
lined above (sections |3.1.H3.1.4[ ) provides a model with no 
consideration of dust effects. This is hereafter referred to 
as model N. Model N provides a basis for comparison with 
photometry corrected for dust attenuation, as described be¬ 
low (sections |3.2.1|3.2.3 1. Panel a) of Fig. shows the g — r 
colour distribution of EAGLE galaxies for this model. 


3.2 Dust model 

We develop a simple empirically calibrated model for dust 
absorption, as opposed to a more physical modelling using 
ray-tracing, which we will present elsewhere. One advantage 
of such a model is that we can easily disentangle the effects 
of dust from those of the SPS modelling on galaxy colours. In 
addition, if we model dust attenuation using galaxy param¬ 
eters that are provided by EAGLE but can also be inferred 
through observation, then we may calibrate the reddening 
of EAGLE galaxy colours empirically to reproduce observed 
trends. Keeping our dust model parametrisation indepen¬ 
dent of certain quantities, such as the gas distribution, also 
allows us more freedom to investigate the extent to which 
certain assumptions affect galaxy colours. 

In our modelling, dust corrections are applied as a mul¬ 
tiplicative factor that reduces a given broad-band luminos¬ 
ity. This factor is estimated at the effective wavelength of 
each filter (SDSS filter parameters taken from 


Doi et al. 


2010, UKIRT hlter parameters taken from |Hewett et al. 
2006), for a given dust prescription (neglecting changes in 
the dust opacity within a broad-band filter is a good approx¬ 
imation provided the dust model has a smooth wavelength 
dependence). In this way, dust obscuration depends on the 
subhalo properties of a galaxy alone, and we avoid handling 
entire SEDs and generating a new interpolation grid for each 
galaxy. This makes the analysis process very efficient. The 
approximation that reddening can be applied after the spec¬ 
trum has already passed through a filter affects the g — r 
optical colours by < 0.02 mag over the whole interpolation 
grid for the constant optical depth model discussed below. 


3.2.1 Model GI: Galaxy-independent dust model 

We begin by discussing the simplest dust model introduced 
by Chariot & Fall ( 2000| , hereafter referred to as CF. This 
model includes two contributions to the dust optical depth 
in a galaxy (fd): (*) a transient component due to dust in 
stellar birth clouds (fbc), and (ii) a constant dust screen 
that represents dust in the ISM (riam). The transmission T 
of this model is 
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exp I - [fbc + Tism] () , for t' ^ idisp , 

n\t') = { ) , ._oV ^ (2) 

exp ( -Tiam f ^ j 1 , for t' > tdisp • 

Here, t' is the stellar particle’s age, tdisp is the dispersal 
time of the stellar birth cloud, A is the wavelength of light, 
and fiam and fbc characterise the total dust optical depth 
(due to ISM and birth-cloud, respectively), at wavelength 
A„. When a model SED of a galaxy is fit to an observed 
galaxy, fiam and fbc can be used as htting parameters to 
model dust (e.g. da Cunha et al.|2008 l. Alternatively, these 
parameters can be assigned constant values to model dust 
for a given population of galaxies as in |Bruzual fc Charlot| 
(20031, but this does then not allow for variations between 


galaxies. 

As a first approximation we simply take fism and fbc to 
be constants. 


'T'hc 

’7”ism 


= 0.67 

= 0.33, 


(3) 


with Ac = 5500 A and tdisp = 10^ yr, which were calibrated 
to ht starburst galaxies and were used for the recent analysis 
of the ILLUSTRIS simula tions ([Vogelsberger et al.||2014[ ) by 
Genel et al. (20141 and Torrey et al. (20141. With the op¬ 


tical depths fixed, the colours of an EAGLE galaxy will only 
depend on the SSP modelling. Such an approximation was 
also used in the development of the galaxev model, where 
it was shown to work well when SED fitting a subset of the 
SDSS survey at z = 0.1 ( [Bruzual Sz Charlot||2003| |. The ef¬ 
fect of this simple dust model on gr — r colours can be seen 
by comparing panels a and b in Fig. and is discussed in 
more detail below. 

In this simple model (model GI), the strength of the 
two screen components are hxed for each galaxy (and hence 
do not depend on e.g. its gas mass or metallicity) and are 
also independent of orientation. We discuss refinements of 
the dust model next. 


3.2.2 Model GD: ISM-dependent dust model 

To account for variations in metal enrichment in the inter¬ 
stellar medium (ISM) of galaxies, we use the mass-weighted 


SPH-smoothed metallicity (Wiersma et al. 2009b I calcu¬ 


lated for star-forming gas in each eagle subhalo, Zsp. This 
metallicity calculation is chosen to imitate observational 
measurement techniques ( Tremonti et al.|2004 Zahid et al. 
2014). As demonstrated in S15, the mass-metallicity rela¬ 


tions in EAGLE are significantly affected by resolution, with 
the Recal-25 simulation showing better agreement with the 


observed Zg — M* relation_Jrom|Tremonti et al.|(|2004|) than 


Ref-100 for M* <C 10^°M(^ Near the knee of the mass func¬ 
tion, however, the EAGLE mass-metallicity relation agrees 
well with observation (see S15). 


® Note, however, that the observed mass-metallicity relation suf¬ 
fers from calibration uncertainties that exceed the difference be- 
tweenj^emonti et al.|||2004l) and Ref-100 ([Kewley &: Ellison 2008[l 


and that the more recent re-analysis by Zahid et al.| | |20T4 1 falls 
in between Ref-100 and |Tremonti et al. ||2004|| (see figure 13 of 
S15). 


We take the analytic expression for the M* — Zg mass- 
metallicity relation at z = 0.1 of Zahid et al. (2014), 
Zzn{Mi,), eval uated at the Milk y Way stellar mass, Mmw = 
6.43 X 10^“ Mo ( |McMillan|20TI| |, as the ISM metallicity rep¬ 
resented by the optical depth values of Eq. §. Assuming 
optical depth is proportional to metallicity, we then scale 
the optical depths fbc and tism that appear in Eq. ([§ by 
the factor 


^bc 


ZsF 


Zzi4(-M* = Mmw) 
.^SF 

Zzn{Mi, = Mmw) 


‘The 


(4) 


for each EAGLE galaxy. 

Making the dust optical depth depend on metallicity is 
physically well motivated, as it is indicative of the fraction 
of ISM mass in dust grains. Therefore, we must also take 
account of the gas mass to quantify how much dust is avail¬ 
able for obscuration. We do so by making the dust optical 
depth dependent on the cold gas mass - but still neglect how 
that gas is distributed. 

We approximate the cold gas mass, Mism, by the mass 
in star-forming gas, which in EAGLE means gas above a given 
metallicity-dependent density threshold and below a given 
temperature (see S15). We then scale the birth cloud and 
ISM dust optical depths by the ratio of Mism for the galaxy 
over the value for the Milky Way (which we take to be 10 per 
cent of M*; Mo et al. 2010 McMillan 2011). This scaling 


is derived from observations indicating that optical depths 
approximately scale with the overall gas surface density (Eg) 
of galaxies (e.g. [Grootes et al.|[^13| |Boquien et al.|[^013p . 
By taking Eg as approximately oc Mism/R*^, and since the 
z = 0.1 stellar mass-size relation is relatively flat for both 
EAGLE (S15) and observed galaxies ( |Shen et al.||2003| ), we 
approximate that r oc Eg oc Mism. Neglecting the 
dependence maintains a relatively simple parametrisation, 
and appears to have little effect on galaxy colours, due to 
the limited mass range over which reddening is signiheant. 

Such a scaling has the desired effect of reddening gas- 
rich spiral galaxies more than gas-poor elliptical galaxies 
at the same M*. The galaxy g — r colour distributions for 
the model including metallicity and gas fraction dependent 
reddening (model GD) are shown in Fig.[^. For comparison, 
we also show the g—r colour distributions for EAGLE galaxies 
where the value of Zzi 4 (M*) is used in Eq. (|^, instead of 
Zsp. Because low Mism values provide low optical depths 
for the stellar mass range where EAGLE and observed mass- 
metallicity relations differ (M* <C lO^^M©), both treatments 
produce similar distributions. 


3.2.3 Model GD-t-O: ISM-dependent dust model with 
orientation effects 

Finally, we represent the dependence of obscuration on ori¬ 
entation with a simple toy model. We assume the dust ge¬ 
ometry to be an oblate spheroid, with major to minor axial 
ratio of g = a/b — 0.2. This q value is commonly used to 
represent an axial ratio typical of the intrinsic stellar distri¬ 
bution in disc galaxies (e.g. Tully fc Fisher|1977 1. We assign 
to each EAGLE galaxy an orientation w = cos{6), where 9 
is the angle between the galaxy’s minor axis and the line of 
sight. To obtain a random orientation we randomly sample 
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w from a uniform distribution over the interval [—1,1]. The 
line-of-sight depth through the dust spheroid is then 


l{w) 


y/q^ + (1 — 


(5) 


We then scale tism as 


^ism 


iM 

(i) 


( 6 ) 


so that the mean optical depth does not change. This scaling 
reduces the amount of dust obscuration for most galaxies by 
a small amount, yet increases r by a factor ~ 2 for a small 
number of ‘edge-on’ systems. We assume a value of q that is 
appropriate for discs, but we note that elliptical galaxies - 
provided they have little cold gas - are not strongly reddened 
anyway, hence this orientation correction is not important 
for them. The g — r colour histograms including orientation 
effects are shown as model GD-I-O in Figure Si. We also 
show the colour distributions produced using more oblate 
geometries, with axial ratio values g = 0.1 and q = 0 . 02 , for 
comparison. 


4 RESULTS 

In this section we examine the effects of dust modelling on 
the colours, luminosities, and colour-magnitude diagrams of 
EAGLE galaxies taken from the Ref-100 and Recal-25 models 
at redshift 2 = 0 . 1 . 


4.1 Galaxy colours as a function of stellar mass 


We contrast g — r colours for EAGLE galaxies in narrow 
(0.3 dex) bins of M* for different models of dust absorp¬ 
tion in Fig. (panels a-d). In all panels, the blue line cor¬ 
responds to the observed distribution from the GAMA sur¬ 
vey ( Taylor et al.|20l4 |. Different panels show models with 
no dust (model N, panel a), a dust model that is indepen¬ 
dent of galaxy type (model GI, panel b), a model in which 
dust opacity depends on metallicity and gas fraction (model 
GD, panel c), and finally a model that in addition accounts 
for orientation effects (model GD-fO, panel d). All mod¬ 
els are shown as histograms, normalised to have unit inte¬ 
gral. Models including dust are plotted as solid histograms, 
while the dashed histograms represent model N in all panels. 
Ref-100 distributions for the fiducial dust models are plotted 
in black, whereas Recal-25 distributions are plotted (for the 
lowest mass bin) in red. Model variations are also plotted for 
Ref-100, with a model using the observed mass-metallicity 
relation in panel c shown in green (see section |3.2.2[ ) and 
models with alternative q values in panel d shown in green 


and yellow (see section 3.2.31. 

The observed g — r colours shift from very red massive 
galaxies (usually termed the ‘red sequence’), to a broader 
distribution in colours for 10 ®'® < M^/Mq < 10 ^®'®, and 
Hnally a blue population (usually termed ‘blue cloud’) at 
lower stellar masses. There is little evidence for a strong ‘bi¬ 
modality’ in observed colours even though the data is often 
interpreted that way. Such an interpretation is perhaps due 
to the bimodality seen at a fixed optical magnitude, where 
blue galaxies are pushed into higher luminosity bins. 


Before comparing the models to the data, we inves¬ 
tigate the effects of dust modelling, going from high- to 
low-mass galaxies (top to bottom rows in panels a-d). The 
most massive galaxies in the EAGLE Ref-100 simulation 
(M*/M 0 > 10 ^^'®, top rows) have a relatively extended 
intrinsic colour distribution (model N). Including a model 
with dust reddening independent of galaxy properties incor¬ 
rectly reddens the reddest galaxies even more (model GI) 
but taking into account the relatively low cold-gas masses 
of these galaxies returns the colours to close to their intrin¬ 
sic values (model GD). The tail of bluer massive galaxies is 
signihcantly affected by dust, yielding a mono-modal distri¬ 
bution in the highest-mass bin. Including orientation effects 
(model GD-fO) gives a slightly broader colour distribution 
for the fiducial axial ratio value of g = 0 . 2 . 

A similar trend is noticeable for galaxies in the sec¬ 
ond most massive bin (10^®'® < M»/Mq < 10^^'®, second 
row from the top). Though the hducial value ot q — 0.2 for 
GD-fO produces a similar distribution to GD in this bin, 
varying q values appears to have the strongest effect here. 
Smaller q values show a more pronounced bimodality, as the 
majority of galaxies are subject to less reddening, while a 
minority of ‘edge-on’ galaxies are heavily reddened to g — r 
colours > 0 . 8 . 

The third bin (10®'® < M*/Mq < 10^®'^, third row 
from the top) also behaves similarly: the intrinsically bluest 
galaxies get reddened slightly more than the intrinsically red 
galaxies with scatter due to orientation having a negligible 
effect. 

Finally, the colour distribution for galaxies in the least 
massive bin ( 10 ®'^ < M*/M 0 < 10 ®, bottom row) is also 
shown for the Recal-25 simulation. There is a large differ¬ 
ence in colours between Ref-100 and Recal-25 for the least 
massive galaxies, which is predominately a resolution ef¬ 
fect: at the resolution of Ref-100 the star formation rates 
in these low-mass galaxies is underestimated (S15) which 
makes the simulated galaxies too red. This striking resolu¬ 
tion dependence is not surprising. In Ref-100, galaxies of 
mass M* ~ 1 O®M 0 are represented by only ~ 10 ® star par¬ 
ticles, and for a typical cold gas fraction of 10 per cent, by 
only 100 star-forming gas particles. We demonstrate in Ap¬ 
pendix that red and blue sequence colours for galaxies 
across the 10 ®'® < M*/M 0 < 10 ®®'® rang^are quite simi¬ 
lar in Ref-100 and Recal-25 - which gives us conhdence that 
Recal-25 gives numerically converged answers for the bottom 
row of Fig. However, the different environments probed 
by the Ref-100 and Recal-25 models also contribute to the 
difference in colours, in particular the strength of the red 
sequence, because the larger volume contains a population 
of satellite galaxies in massive halos. This is also shown in 
Appendix and discussed further in section Taking into 
account dust obscuration and orientation effects has little 
effect on the colours in Recal-25 for these low-mass galaxies, 
with model GD+0 and N yielding nearly identical colour 
distributions. 

We now turn to comparing the colours of EAGLE galax¬ 
ies to the data, going from top (most massive) to bottom 
(least massive) bins in stellar mass and focusing on model 


A larger mass range is used in Appendix ^ than in the third 
row of Fig.niso that Ref-25 and Recal-25 are sufficiently sampled. 
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a) No dust (N) b) Galaxy-iiidepeiident model (GI) 




c) Galaxy-dependent model (GD) d) Galaxy-dependent model with orientation effects (GD-l-0) 




Figure 1. Rest-frame g — r colour distributions for EAGLE galaxies at redshift 2 : = 0.1, using 4 different models {panels a-d) for four non¬ 
contiguous ranges in stellar mass as indicated in the legend {top to bottom). Black lines indicate the fiducial Ref-100 galaxy population 
while red lines indicate the higher-resolution Recal-25 simulation. Dashed lines denote the unobscured SED (model N); these are repeated 
in each panel for comparison. Models shown are: model N without dust (panel a, see ^3.1.5[ |) model GI with galaxy independent dust 
(panel b, see § |3.2.1| l, model GD where the dust obscuration depends on gas fraction and metallicity (panel c, § |3.2.2[ |, and model GD-hO 
that in addition takes into account orientation effects (panel d, § |3.2.3[ |. Green and yellow lines show model variations in panels c and d 
(see sections |3.2.2[ 3.2.3| for details). Blue lines represent observed galaxy colours for the volume-limited sample of GAMA galaxies from 
|Taylor et al.| l |2014 l. The figure shows the subtle quantitative effects that our different dust models have on the colour distributions of 
EAGLE galaxies in various M* regimes, as is discussed further in the text. 
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Normalised Weighted Point Density 
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Figure 2. Rest-frame g — r colour-stellar mass diagrams for EAGLE galaxies using photometry and dust reddening from model GD-l-0 
at 2 = 0.1 for simulation Ref-100 (left panel), and a composite set of EAGLE galaxies obtained by combining the higher-resolution model 
Recal-25 for galaxies with M* < 10^ Mq and Ref-100 for M* > 10® Mq {right panel), cross-fading the sampling probability of the 
two galaxy populations linearly in \ogiQ{M*/M q) (see text for details). Colours represent the point density of EAGLE galaxies of given 
M* and g — r colour (see text for details), with black points representing individual outlying EAGLE galaxies. Contours represent the 
colour-M* distribution for a volume-limited set of GAMA galaxies from [Taylor et al.| | |2014||, with gre y points representing individual 
outlying galaxies. The masses of observed galaxies are obtained through SED fitting, see [Taylor et al!] l |2014| |. The colours are mapped 
to their equivalent contours in the colour bar. The colour bar covers 2 dex in point density with a contour spacing of 0.28 dex. A lack of 
numerical resolution makes lower-mass EAGLE galaxies too red in the left panel. The transition of colours of observed galaxies, from the 
red-sequence at stellar masses above M* ~ to the blue cloud at lower stellar masses, is reproduced in the simulation, although 

the blue cloud extends to slightly higher M*. 


GD+0, Fig. Bi)- At the massive end, the observed red se¬ 
quence galaxies are about 0.05 mag redder in the data than 
in EAGLE. As the optical colours of old (> 10 Gyr) stel¬ 


lar populations are dominated by metallicity effects | Bell fc 
Rodgers|1969 this small colour difference is attributable to 


SSP metallicities. The AL* — relation for EAGLE galaxies 
(S15) is seen to lie slightly below (by less than 0.1 dex) ob¬ 
servational data in galaxies with Af*/M 0 > 10 ^^, resulting 
in a slightly bluer red sequence colour. The data also has a 
tail to even redder colours not present in EAGLE. In contrast, 
the most massive EAGLE galaxies have a tail to bluer colours 
resulting from recent star formation. It could be that such 
star formation is shielded more effectively in the data (ie. 
the value of fbc used is too low), or alternatively that our 
AGN feedback scheme does not quite suppress star forma¬ 
tion sufficiently. The higher than observed gas fractions for 
galaxy clusters in EAGLE (S15) could also contribute to the 
enhanced SFR of some simulated BCGs. 


The red sequence of galaxies with < AL*/Mq < 

10 ^°'® is very similar in the data and the simulation, but 
in EAGLE there are significantly more blue galaxies. The 
blue cloud starts to appear in the data for galaxies with 
10 ®'® < M*/M 0 < 10 ^®'^, and its colour is very similar in 
EAGLE. However, in EAGLE the blue peak is stronger and the 
red peak occurs at a slightly redder colours {g — r = 0.75 
compared to the observed value of 0.7). Using smaller values 
of q in GD-I-O does not improve agreement with observation 


here. Dust reddening and orientation effects already play lit¬ 
tle role in setting the colours of EAGLE galaxies in this mass 
bin. Finally, in the lowest mass bin, 10®'^ < AL*/M 0 < 10®, 
there is excellent agreement in the colour distributions of 
simulation and data and once more our dust reddening mod¬ 
els are unimportant in setting EAGLE colours. 

This level of agreement between galaxy colours in the 
simulation and the data is encouraging. By including metal¬ 
licity and orientation effects in our dust treatment, we pre¬ 
vent the significant colour shift seen in the simple GI model. 
The validity of our dust model is discussed further in sec¬ 
tion Despite the good agreement, there are some clear 
discrepancies between eagle and the GAMA colour distribu¬ 
tions. These can be seen in the widths and relative strengths 
of red and blue populations. The latter discrepancy reflects 
the finding of S15 that the transition from actively star¬ 
forming to passive galaxies occurs at slightly (by a factor of 
~ 2 ) too high mass in eagle. 

The dependence of galaxy colours on stellar mass is 
further illustrated in Fig. where the number density of 
galaxies in EAGLE with given rest-frame g — r colour (com¬ 
puted using model GD-I-O) and stellar mass is compared to 
a volume-limited sample of GAMA galaxies taken from|Taylor| 


et al. (2014). Results are plotted down to stellar masses of 


10 ®''M 0 , below which volume corrections due to the influ¬ 
ence of line-of-sight structure become increasingly uncertain 
in the data (Taylor et al. 20141. The colour bar shows the 
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point density of EAGLE galaxies and how these map to the 


Taylor et al. (20141 contours. 


The left panel shows the galaxies taken from simula¬ 
tion Ref-100. The simulation reproduces the trend seen in 
the data from galaxies being red above a stellar mass of 
M* ~ 1 O^° ®M 0 to being predominantly blue below that. 
However, as also seen in the previous figure, there is a pop¬ 
ulation of red (<; — r ~ 0.7) low-mass {M^, ~ IO^Mq) galax¬ 
ies in EAGLE that is not seen in the data. These galaxies 
are modelled using only ~ 1000 star particles; the right 
panel of Fig. [^therefore uses the higher resolution simula¬ 
tion Recal-25 for galaxies below IO^Mq and Ref-100 above 
IO^Mq, cross-fading one simnlation into the other linearly 
in log(M*f] With this we aim to show the colour-mass dis¬ 
tribution for a larger range of well-resolved galaxies (> 1000 
star particles), while avoiding a discontinuity that renders 
the overall distribution less clear. 

Combining these two resolntions, the colours of EAGLE 
galaxies at given M* track the data from the GAMA galax¬ 
ies ( Taylor et al.|20l4 l well. Both display a red sequence of 
massive galaxies which becomes redder with increasing stel¬ 
lar mass, and g — r ~ 0.7 at M^, ~ 10^°'®Mq. The simnla¬ 
tion also reproduces the width of that sequence, albeit with 
a shallower slope. A bine cloud of galaxies appears both in 
EAGLE and GAMA below M* ~ IO^^'^Mq, with p — r ~ 0.45 
at Mi, = 1 O^°'®M 0 . At decreasing stellar mass, the loca¬ 
tion of the blue cloud becomes bluer, reaching g — r ~ 0.35 
at Mi, = 10® M 0 . Overall we find that EAGLE reproduces 
the mean trends in galaxy colours well. Though the erad¬ 
ication of the faint red sequence in this sample is at least 
partly due to improved sampling, it also comes about be¬ 
cause these galaxies are much less abundant in the higher 
resolution Recal-25 simulation (further discussion of the ori¬ 
gin of this faint red population can be found in Sectionj^and 
Appendix E- 

In addition to the mean location of galaxy colours, there 
are outliers in both data and simulation. The GAMA data 
display a scatter to extremely red colours {g — r > 1 ) at 
all stellar masses only seen for one high-mass outlier in EA¬ 
GLE. There is also a scattering of galaxies ~ 0.1 mag bluer 
than the main locus in GAMA that appear in EAGLE as well. 
Finally, EAGLE has some very massive, relatively blue galax¬ 
ies {Mi, ~ 1 O^^ ®M 0 , g — r ^ 0 . 6 ); although there are such 
galaxies in GAMA as well, they are more numerous in EAGLE, 
as is more easily seen in Fig. We suggested before that 
these either imply too little dust reddening in star forming 
regions in EAGLE, or simply that some of these massive EA¬ 
GLE galaxies are undergoing too much star formation despite 
the inclusion of AGN feedback. 


® The sampling frequency of Recal-25 galaxies in this plot are 
weighted a factor of 64 higher than Ref-100 galaxies to account for 
the smaller volume and a further factor of 2 lower to account for 
the boosted number counts in Recal-25 caused by poor sampling 
of large scale power in the smaller volume. This means that for 
masses < 1 O®M 0 each Recal-25 galaxy contributes a factor of 32 
more to the point densities in Fig. [^than a Ref-lOO galaxy at 
mass > 1 O®M 0 . 


4.2 Luminosity functions and colour-magnitude 
diagram 

Luminosity functions for model GD+0 in rest-frame 
ugrizYJHK broad-band filters are plotted using absolute 
AB magnitudes in Fig. The simulations Ref-lOO and 
Recal-25 at redshift z = 0.1 are shown with Poisson error 
bars as solid black and red histograms respectively, becom¬ 
ing dashed when there are fewer than 10 galaxies per ~ 0.6 
mag bin. For Recal-25 the bins are correlated, as can be seen 
for example in the u-band for bins M — 51og(h) = —17.5 
to —18.5, due to poor sampling of large-scale modes in the 
small volume. There is generally good agreement between 
the two runs, with Recal-25 typically less than a factor of 
two (0.3 dex) higher at the faint end, and by mnch less for 
the redder bands. Note that this higher-resolution simnla¬ 
tion does not sample the exponential cut-off at high lumi¬ 
nosities because of its small volume. Differences in resolution 
are most noticeable in bluer colours, particularly in u. As dis¬ 
cussed before, in small galaxies the stellar feedback events 
driving outflows are poorly sampled and the star forming 
components are poorly resolved. As a result star formation 
rates and thus intrinsic colours are subject to considerable 
resolution effects. We see that the higher-resolution simu¬ 
lation yields higher star formation rates and hence bluer 
colours. We again note that this is not just a resolution is¬ 
sue: the Ref-lOO volume contains a population of faint red 
quenched satellites of massive galaxies, which are simply not 
present in the much smaller Recal-25 volume, as discussed 
further in Appendix]^ LFs in longer wavelength bands are 
consistent in shape between the two simulations (despite the 
small volume simulation being noisier). 

Observed luminosity functions from |Loveday et al.| 


(20121 are plotted in each of the ugriz bands, which we 


fit with a single Schechter (1976 1 function. 


1 


dn 

dL 


exp(-L/L*) 


dL 


(7) 


Single Schechter function fits are also taken from |Driver| 
( |2012[ ) and [Jones et al.| ( |2006| ). For the ugriz bands these 
are shown as grey shaded regions which are bounded by the 
fits to the observed luminosity function of [Loveday et al.j 
( (M2| ) and|D riverj ( |2012| , both based on data from GAMA. 
The differences between these observed luminosity functions 
result from the use of Kron and Petrosian magnitudes, re¬ 
spectively. The thickness of the grey band is thus a measure 
of how these different aperture choices affect the Schechter 
fit. For the YJHK band, we plot published Schechter fits, 
which are based on ukidss data. 

There appears to be some discrepancy between 
Schechter fits from the observational papers and directly ob¬ 
served luminosity functions in the optical, as can be seen by 


comparing the Loveday et al. (20121 data with the Schechter 


fits. In particular, the data points appear systematically 
higher than the Schechter fit at the faint end and below 
the fit at the bright end. This is most visible in the «-band 
where the shaded region is narrowest. This could be a conse¬ 
quence of intermediate magnitude bins dominating the fit as 
this is where observational errors are minimal. It also shows 
that the single Schechter function is not a good fit to the ob¬ 
servational data. In particular, it is unable to represent the 
observations at the faint end accnrately (e.g. Loveday|1998|). 
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M - Slogj^g h [mag] 
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Figure 3. Luminosity functions in each band of the ugrizYHJK photometric system plotted for the fiducial 100 Mpc simulation 
(Ref-100) in black and the re-calibrated high resolution 25 Mpc simulation (Recal-25) in red, both at redshift 2 = 0.1, using the GD-l-0 
dust model. The Ref-100 function is plotted down to the faintest magnitude bin at which most galaxies are represented by >100 star 
particles. The simulation results from ugriz SDSS bands are compared to the Petrosian GAMA survey luminosity functions taken from 
Loveday et al. /c-corrected to rest-frame magnitudes and plotted as empty circles. The region bound by the Schechter fits to the 

Driver] < |2012| | and [Loveday et al.| l |2012| | luminosity functions is shaded in grey. As these two luminosity functions are measured using 
Kron and Petrosian magnitudes respectively, the grey area indicates the difference due to aperture definition. Schechter fits for UKIDSS 
band luminosity functions are taken from GAMA (Driver et. al. 2012) and from 6DF-|-2MASS < | Jones et al.|[200^ where available. All 
magnitudes are in the AB system. Error bars refiect Poisson errors with dashed lines indicating bins containing < 10 simulation galaxies. 
The EAGLE luminosity functions are similar to the observed fits across the spectral range, with some discrepancies discussed in Section]^ 
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Model 

Description 


0* [h^ cMpc 

r-band 

a 

M* - 51ogig(/i) [mag] 

N 

GI 

GD 

GD+O 

No dust 

Galaxy-independent dust model 

Z and Mism dependent 

GD with orientation dependence 

T.olp g X 10-3 

7.6lo g X 10-3 
8.2t3 0 X 10-3 

S.Slg 1 X 10-3 

-I qq-I- 0.04 
-‘-•^'^-0.02 
-| Qq-|-0.05 
J--^'5 _o.02 

1 Qq-|-0.04 

02 

01 r\-\-0.‘2 
^-*^•^-0.2 

-20.7to\ 

Data 

Loveday et al. 

(2012| 

9±0.7 X 10-3 

-1.26±0.01 

-20.7±0.03 


Table 2. Best-fitting Schechter function (Eq.[^ parameters for EAGLE AB-magnitude luminosity functions in the r-band for different 
dust models in simulation Ref-100 at redshift 2 = 0.1, and the observed luminosity function from [Loveday et al.| l [^12| |. The EAGLE 
and observed luminosity functions are all fit over the magnitude range —23.2 < r < —14.2, Icr errors on the best-fit parameters were 
computed using jackknife sampling. 


The Driver (20121 and Jones et al. (20061 fits agree reason¬ 


ably well, except for the faint-end slope in the Tt'-band. The 
single Schechter fits are used here as simple indicators of 
the shape, position and normalisation of the observed lumi¬ 
nosity functions, but clearly their exact location depends on 
details of how galaxies are identified in the data, and pos¬ 
sibly on the range and assumed errors used in the fitting 
procedure. We compare the parameters of Schechter fits to 
EAGLE luminosity functions to observational fits in Table 

From Table [ 2 ] we see that the dust treatment has little 
effect on the shape of the r-band luminosity function. The 
effect of including dust using the GI model makes the knee 
position, M*, 0.5 mag fainter magnitudes and increases the 
normalisation, cf),. Scaling dust absorption by galaxy prop¬ 
erties in GD serves to increase </)* further, and makes M* 
brighter by 0.3 mag. The GD and GD-I-O model luminosity 
functions provide Schechter fits that agree with the obser¬ 
vational values within the errors. The faint-end slope, a, 
shows very little variation between all dust models, and all 
are consistent with the observational fit within the errors. 

Comparing the GD-I-O EAGLE luminosity functions to 
the data in Fig. shows a striking overall consistency from 
the UV to the NIR bands. The deviations are mostly of the 
same order as differences in fits to the published luminosity 
functions of different authors. The agreement is particularly 
good in the optical bands ugriz, where EAGLE tends to fall 
mostly inside the grey band that represents the dependence 
of the luminosity function on the choice of aperture. The 
excellent agreement over such a wide range of colours sug¬ 
gests that EAGLE forms the correct number of galaxies of a 
given stellar mass and that those galaxies have realistic star 
formation histories and metallicities. 

Comparing blue bands (u—g) to redder bands (J—K) at 
the faint end, we notice that the EAGLE Ref-100 luminosity 
functions tend to be slightly low in blue bands relative to 
the data, but high in the red bands. This is a consequence 
of Ref-100 producing slightly too many low-mass galaxies 


(S15) which have too low star formation rates (Furlong et al. 


2014). Resolution also plays a role: we plot galaxies with 


more than 100 star particles, where we know that the stellar 
feedback events generating outflows and star formation rates 
at the faint end are poorly resolved. In addition, even fainter 
galaxies with high star formation rates cannot scatter into 
the faint-end bins since we impose a cut in mass and not in 
magnitude. 

The EAGLE luminosity function tends to drop below the 
observations at the ‘knee’ (L*) in the Schechter function. 


particularly in the bands red-ward of r. This is consistent 
with a slight underestimate in the masses of more massive 
EAGLE galaxies, as seen in the mass function plotted in S15. 

The JHK bands also appear to have generally some¬ 
what steeper faint-end slopes (parameter a in Eq 0 in EA¬ 
GLE than the Schechter fits to the data. However, the data 
itself also shows a upturn at the faint end relative to the ob¬ 
served Schechter fits (circles in Fig.j^ see also Driver|2012 1 . 
Generally, the single Schechter function fit tends to under¬ 


estimate the luminosity function at the faint end (Loveday 


1998 1 . This is more pronounced in the JHK bands where the 


NIR sky is relatively bright (e.g. Sivanandam et al. 2012), 
leading to large uncertainties in the faint end data. 


4.3 The g — r colour-magnitude diagram 

The colour-magnitude diagram of model GD-I-O for EAGLE is 
plotted in Fig.[^ As in the right panel of Fig. we combined 
faint galaxies from the higher-resolution simulation Recal-25 
for galaxies with M* < 1 O®M 0 with galaxies from simulation 
Ref-100 at higher masses. As before, colours represent the 
number density of EAGLE galaxies in this plane, whereas con¬ 
tours show the corresponding data from the volume-limited 


catalogue of GAMA taken from Taylor et al. (20141. For both 


simulation and data, we only show galaxies with stellar mass 
M* > 1O®'^M0. 

Fig. |4] now contains only ‘observable’ quantities for 
GAMA galaxieB and in particular does not require any SED 
fitting. The overall agreement between EAGLE and the data 
is very good: the location and slope of the red-sequence 
is in good agreement, the transition to the blue cloud at 
Mr — 51og(/i) ~ 20 is reproduced, and the blue g — r ~ 0.3 
colours of the fainter blue galaxies are also reproduced. At 
the imposed mass cut of M* > 1 O®’^M 0 , the colours and 
magnitudes of EAGLE galaxies agree well with the data. The 
bright blue galaxies that appear in EAGLE become slightly 
more discrepant with the data in this plot than in Fig. 
due to the fact that bluer galaxies generally possess brighter 
r-band magnitudes for the same M*. 

The level of agreement between simulation and data in 
the colour-magnitude diagram shown in Fig. is remark¬ 
able and suggests that these EAGLE simulations provide a 


° In practice the diagram still depends to a small extent on the 
applied cut in stellar mass at the faint end and on the choice of 
aperture to measure magnitudes. 
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Figure 4. Rest-frame g—r colour as a function of r-band absolute 
magnitude for EAGLE galaxies (colours) compared to a volume- 
limited sample of GAMA galaxies i jTaylor et al.| [20141 contour 
lines), contours and point shading is the same as in Fig EA¬ 
GLE photometry is obtained using the GD-hO model. A compos¬ 
ite EAGLE galaxy population is used, consisting of galaxies from 
Ref-100 at M* > IO^Mq and Recal-25 at M* < IO^Mq as in 
the right panel of Fig[^ There is excellent agreement between 
EAGLE and GAMA in the location and slope of the red-sequence 
(g — r 0.7 at Mr — 5 log^Q h = —20), the appearance of a blue 
cloud of galaxies with g — r 0.45 at that magnitude, which be¬ 
comes increasingly blue {g — r 0.3) for the fainter galaxies with 
Mr — 5 log h —17.5. 
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relatively realistic population of galaxies at low redshift, 
and that the modelling of emission and dust obscuration 
in model GD+0 works well. 


5 DISCUSSION 


The EAGLE simulations were calibrated to reproduce the lo¬ 
cal GSMF and galaxy sizes by appropriate choice of the pa¬ 
rameters in the subgrid model for feedback. As stellar mass is 
closely linked to NIR luminosity, the consistency of GSMF 
proxies such as the A-band luminosity function shown in 
Fig. [I is not surprising, at least at 2 ; ~ 0. However, con¬ 
sistency with the luminosity functions in other broad band 
filters is not automatic because stellar mass, star formation 
history, metallicity and dust obscuration all play a role. 

In Section we focused on our GD and GD+0 pho¬ 
tometric models where dust absorption is approximated by 
a simple two component screen, with optical depths that 
vary with galactic gas content, metallicity and orientation 
(GD-l-0 only). The colour distributions as a function of stel¬ 
lar mass and r-band magnitude for model GD-I-O in Figs.[^ 
[^andj^show a level of agreement between simulated galaxies 
and observations that appears unprecedented for hydrody- 
namical simulations, and comparable to that achieved for 
semi-analytical models (e.g. [Gonzalez et al.|2009{|Henriques| 
et al.|2()14 |. The EAGLE luminosity functions also agree well 


with observations over a range in wavelengths from optical 
to NIR (Fig. 1^. The relatively good agreement for number 
density, luminosity and colour, suggests that in EAGLE each 
dark matter halo forms a galaxy with stellar mass, age and 
metallicity close to those inferred from observation. The sim¬ 
ilar level of agreement in observed colour-magnitude space 
also rules out a potential circularity resulting from using the 
same photometric model to infer stellar mass from observa¬ 
tional data as is used in the simulations. Even though the 
overall level of agreement is good between EAGLE and the 
data, there are discrepancies. 

There is an excess of bright (Mr — 5 logjQ h < —20) blue 
(g — r < 0.6) EAGLE galaxies relative to data, apparent in e.g. 
Fig-i Such an excess is seen in all our photometric models 
(see Fig. [^, but is least apparent for model GD-fO where 
the recent star formation, that is the root cause of the blue 
colours, is most strongly obscured by dust. This parallels the 
findings of S15, that the fraction of passive EAGLE galaxies 
is too low at the high-mass end, relative to observations. 
It may be that massive EAGLE galaxies are too highly star 
forming, perhaps as a consequence of insufficient suppression 
of star formation by AGN. Whether the good agreement in 
the colour of the blue cloud, despite the underestimate of the 
median star formation rates at these masses (S15), suggests 
underestimated reddening is discussed further below. 

The EAGLE red sequence, at g — r ~ 0.75, is flatter than 
observed, both when plotted as a function of stellar mass 
(Fig.[2| and when plotted as a function of absolute magni¬ 
tude (Fig. 1^. The flatter slope may be attributable to the 
dependence of stellar metallicity on galaxy mass, Z*(M*); 
although stellar metallicities of EAGLE galaxies agree well 
with the data at the massive end, they fall less rapidly with 
decreasing stellar mass compared to the observational data 
of Gallazzi et al. (20051, as shown in S15. Numerical res¬ 
olution may play a role here, because the Z*(M*) of the 
higher-resolution simulation Recal-25 does agree with the 
data; see S15. 

There is an abundant population of red (g — r ~ 0.7), 
low-mass (M* ~ lO^M©) galaxies in simulation Ref-100 that 
is not observed (see Fig.[^). The comparison of simulations 
Ref-25 and Recal-25 in Appendix B shows that this is at 
least partially due to a lack of numerical resolution. Indeed, 
star formation and outflows driven by feedback in these 
galaxies are poorly resolved and poorly sampled, leading to 
too low values of the specific star formation rate (and corre¬ 
spondingly too high passive fractions) and too high metallic¬ 
ities (see S15). Appendix [A| discusses the re-sampling tech¬ 
nique we use in all models to try to mitigate poor sampling. 
Although this goes some way towards improving the mod¬ 
elling, it does not eliminate the discrepancy. Because the 
re-sampling is a post-processing step, it cannot help with 
the poor sampling of stellar feedback in these low-mass sys¬ 
tems within the simulation. Related resolution problems are 
more intractable, and higher-resolution simulations are re¬ 
quired to alleviate them. 

Comparing simulations Ref-100 and Ref-25 that have 
identical numerical resolution (and the associated poor sam¬ 
pling of star formation in M* ~ IO^Mq galaxies), yet differ 
in simulated volume size, allows us to isolate the effects of 
environment (see Appendix 0. Although on average the 
colours of galaxies agree well between these simulations, the 
presence of faint red galaxies is much more pronounced in 
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the larger volume. This is because many of these galaxies 
are satellites of more massive systems that are absent in the 
smaller volume. The fraction of satellites increases at lower 
stellar masses, and in the range < M* < IO^Mq 

comprises ~ 46% of the galaxy population in the Ref-100 
simulation and ~ 33% in Ref-25. Evidently, satellites con¬ 
tribute significantly to the colour distribution at low masses. 
At present we cannot verify whether improved resolution will 
also reduce the suppression of star formation in, or decrease 
the metallicities of, small satellite galaxies, which would im¬ 
prove the colours of M* ~ 10®Mq galaxies compared to 
data. We conclude that the redder colours of low-mass galax¬ 
ies in EAGLE relative to data is at least partially a result of 
resolution, stemming from poor sampling of star formation 
and feedback. The improved agreement with the data that 
comes about from using a composite sample of Recal-25 and 
Ref-100 galaxies relative to using the Ref-100 sample alone 
is thus mainly due to improved numerical resolution, but 
also to the exclusion of red satellite galaxies that are not 
present in the smaller volume. 

The level of agreement between eagle colours and 
the data also depends on the realism of our dust redden¬ 
ing model. Fig. illustrates how dust reddening depends 
on the assumptions made in models N (no reddening) to 
model GD-I-O (gas metallicity, gas mass, and orientation- 
dependent reddening). Differences between these models are 
typically of order A{g — r) ~ 0.1. A dust model that is in¬ 
dependent of galaxy properties (such as GI) incorrectly red¬ 
dens red galaxies. A reddening model that takes into account 
the gas mass (GD) resolves this inconsistency, with most 
of the remaining effects of reddening affecting blue bright 
galaxies. Overall, we find that the details of the dust treat¬ 
ments make relatively little difference to galaxy colours so 
that differences with observations are more likely due to the 
ages and metallicities of the stars rather than dust obscura¬ 
tion. This may be due to the relatively small effect of dust 
at redshifts 2 ~ 0.1. We leave the investigation of evolution 
of colours and luminosities in EAGLE to a future work. 

Taken at face value, the specific star formation rates 
of star-forming galaxies in EAGLE are lower than inferred 
from observations by ~ 0.2 dex (S15). The fact that the 
colours of those same galaxies nevertheless agree with the 
data might imply that we underestimate dust reddening. 
Indeed, an underestimate of the gas fraction would lead to 
both an underestimate of the specihc star formation rate 


and, at fixed metallicity, the attenuation. Lagos et al. (20151 
have shown that at M* ~ IO^^Mq the median H 2 fraction in 
EAGLE is about 0.2 dex lower than observed, and that this 
discrepancy goes away at higher masses and for the higher- 
resolution Recal-25 model. 

Systematically lower attenuation for faint galaxies could 
thus be attributable to their low gas masses in the simula¬ 
tions. However, more complex models yield non-zero levels 
of attenuation even for very low gas surface densities (e.g. 
Boquien et al. 2013 L The realism of mixed screen models 


as used here has been shown to break down when screens 


are optically thick (e.g. Disney et al. 19891. As dust op¬ 
tical depths are expected to be higher in blue bands, our 
dust prescription may be too crude to reproduce the data 
at higher levels of obscuration. This could contribute to the 
bluer colours of massive galaxies in EAGLE. 

However, it is also possible that the levels of obscura¬ 


tion are realistic, but that star formation rates are overes¬ 
timated in the data due to the absolute calibration of ob¬ 
served tracers of star formation. The calibration of star for¬ 
mation rates from tracers rely on assumptions about the in¬ 
trinsic UV continuum (from population synthesis modelling) 
and absorption at short wavelengths, as well as an assumed 
form for the IMF (e.g. Kennicutt] 1998). The cumulative 
build-up of stellar mass in EAGLE is lower than observed by 
about 0.1 dex, whereas the star formation rate is lower than 


that observed by 0.2-0.4 dex, depending on redshift (Furlong 


et al. 20141. This slight tension may suggest a small overesti¬ 


mate of the observationally inferred star formation rates. Es¬ 
timating intrinsic properties from observables of simulated 
galaxies (such as star formation rates) may help to clarify 


these issues, see e.g. the recent study by Torrey et al. (20141. 

The dust model we developed here was designed to be 
as simple as possible, yet to avoid unrealistic levels of red¬ 
dening. The model assigns a single value of reddening per 
galaxy without taking into account the non-uniform distri¬ 
butions of dust apart from that assigned to birth clouds. It is 
possible to make much more detailed estimates of reddening 


using 3D radiative transfer (RT) calculations (e.g. Baes et al. 
2005 Jonsson et al. ]200^ l. We postpone comparisons of the 


current simple model to those obtained with the radiative 
transfer code skirt (Baes et al.|2005i to future work. 


6 SUMMARY & CONCLUSIONS 

We have calculated broad band luminosities of simulated 
galaxies taken from the EAGLE (Evolution and Assembly of 
GaLaxies and their Environments) suite of hydrodynamic 
simulations (S15, Grain et al. 2015), and compared them 
to observations of the redshift 2 ~ 0.1 galaxy population. 
The model uses simple stellar population modelling based 
on the GALAXEV population synthesis models of |Bruzual fc| 
Chariot (2003). To reduce sampling noise arising from sin¬ 
gle young star particle in poorly resolved galaxies, we use 
the re-sampling procedure described in Appendix]^ In all 
models, galaxy luminosities are found by summing the par¬ 
ticle luminosities within a 30 pkpc radius spherical aperture 


for consistency with previous analysis (S15, Furlong et al. 


2014), which has been shown to mimic Petrosian apertures. 


Absolute magnitudes are presented in the AB system. 

We compare and contrast three models of dust obscura¬ 
tion and to model N which neglects dust. Model GI, inspired 
by Chariot & Fall ( 2000| , includes contributions to the dust 
optical depth from the birth clouds of young stars and from 
a constant dust screen, with parameters that are indepen¬ 
dent of the galaxy properties. Applying a single diffuse dust 
correction to all galaxies incorrectly reddens ellipticals and 
we avoid this with model GD in which dust reddening de¬ 
pends on gas phase metallicity as well as gas mass. Finally, 
model GD-I-O uses a simple geometric model to account for 
orientation effects, which are however small. 

These simple models allow us to investigate the depen¬ 
dence of galaxy colours on stellar metallicities and ages, gas 
metallicities, and dust obscuration. Our main conclusions 
are as follows: 

• The GI dust prescription which applies a reddening that 
is independent of galaxy properties, and was used by e.g. 


Torrey et al. (20141, excessively reddens the red-sequence 
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population of galaxies. As a consequence, g — r colours of 
massive (M* > IO^^'^Mq) eagle galaxies are ~ 0.1 mag red¬ 
der than observed, in spite of having ages and stellar metal- 
licities that are similar ro those inferred(Fig [^). Scaling 
dust optical depths with cold gas mass and gas metallicity, 
as in model GD, is more realistic and improves agreement 
with observation (Fig[^ and[^). 

• The red sequence in EAGLE is ~ 0.1 mag bluer in g — 
r than observed for M* > 10 ^^'^Mq, and has a shallower 
dependence on stellar mass than observed (Fig. [^. This is 
most likely a consequence of the dependence of colour on 
stellar metallicities. 

• The appearance of a faint red sequence in the Ref-100 
simulation run (M* < IO^'^^Mq, 0.6 < g — r < 0.8, see left 
panel of Fig that is not observed, is largely an effect of 
numerical resolution. Star formation and outflows are not 
well resolved in galaxies of such low mass. 

• A ‘blue cloud’ of star forming galaxies appears in EAGLE 
below M* ~ 10^° ®Mq, with g — r colour in agreement with 


the GAMA data from Taylor et al. (2014 1 (left panel of Fig[^. 

• There is an excess of bright {Mr — 5 logj^Q h < —20) blue 

(5 “ ^ 0-6) galaxies in EAGLE relative to the data. This 

may be caused by an underestimate of the reddening in star 
forming regions, or an overestimate of the star formation 
rates in these massive galaxies due to insufficient suppression 
of star formation by AGN. 

• The z = 0.1 galaxies taken from EAGLE transition from 


mostly red {g — r ~ 0.7) above M* 
blue (ranging from ^ — r ~ 0.5 at M* 


Iqiu.&Mq to mostly 
/ 10^° ®Mq becoming 
bluer with decreasing mass to <? —r ~ 0.35 at M* ~ IO^Mq) 
at lower masses, follows the colours of GAMA galaxies from 
Taylor et al. (20141 (see Fig[^. However the blue cloud per¬ 
sists to higher than observed stellar masses, consistent with 
a similar trend in passive fractions shown in S15. 

• The 2 = 0.1 galaxy luminosity functions constructed 
from the EAGLE population agree well with data from UV 
to NIR bands, with differences of the order of the difference 
between using Kron and Petrosian magnitudes in the data 
(Fig.|3f. This level of agreement is similar to the agreement 
between the eagle and observed stellar mass functions. In 
particular, there is a slight underestimate in the number 
density of galaxies close to the knee of the Schechter fit, 
and the faint-end tends to be slightly steeper than observed 
in most bands. We note, however, that the faint-end of the 
luminosity function is uncertain, especially in NIR bands, 
and single Schechter fits tend to underestimate the faint-end 
slope ( Loveday|1998 (. Good agreement was not surprising in 
the NIR where luminosities are dominated by stellar mass, 
whereas the good agreement in other bands suggests that 
the star formation histories and metal enrichment in EAGLE 
galaxies are relatively realistic. 

• The z = 0.1 g — r colour versus Mr magnitude diagram 
for galaxies with M* > IO^Mq yields a level of agreement 
with data that is comparable to that of current semi-analytic 
models (Fig. [4| [Gonzalez et al.|2009 Henriques et al.|2014 |. 
The similar colour distributions of N and GD-I-O photome¬ 
tries (Fig. 01) suggests that the dust model plays only a 
minor role in this agreement. This further attests to the rel¬ 
atively realistic evolution of the EAGLE galaxy population. 


The general agreement in the colour and luminosity of 
EAGLE galaxies and observed galaxies suggests that the sim¬ 


ulated galaxies have similar star formation histories, metal 
enrichment processes, and current star formation rates as 
observed galaxies. This makes the EAGLE suite well-suited 
to investigate the physical processes that shape galaxies 
through cosmic time. 
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APPENDIX A: RE-SAMPLING 


Young stellar populations are optically much brighter and 
bluer than older populations, due to the presence of O, B 
and A type stars. As such, the convergence of optical magni¬ 
tudes and colours depend on wavelength and star formation 
in a complex way. Coarse sampling of young star particles is 
expected to introduce Poisson scatter in the colour distribu¬ 
tion of low-mass galaxies. The standard resolution particle 
mass of nig = 2 x 1 O®M 0 is ~2 orders of magnitude larger 


than individual Hll regions (Zaragoza-Cardiel et al. 2014 
[Relano fc Kennicutt|200^ , so this scatter may be artificially 
high in the simulations compared to observed galaxies. 

To mitigate this effect, we re-sample recent star forma¬ 
tion in the simulation outputs at finer mass resolution than 
the simulation gas particle mass. The re-sampling procedure 
assumes that the galactic star formation rate is constant 
over the previous 100 Myr of galaxy evolution. The general 
procedure is outlined below. 


Gas particles with non-zero star formation rates and 
star particles that formed within the past 100 Myr are first 
identihed. Star formation rates of progenitor gas particles 
are then obtained for the selected star particles. We calculate 
the star formation rates for star particles using the stored 


gas density at birth and initial particle masses (see Schaye 


20041, and assume that the progenitor particles lie on the 


equation of state (see S15). We then take the sum of the star 
formation rates found for the gas and young stars to be the 
galactic star formation rate. 

We then use some discretisation mass, rudis, to repre¬ 
sent the mass of re-sampled star particles. Each selected star 
or gas particle is split into the integer number of equal-mass 
subparticles yielding a subparticle mass closest to mdis. We 
also obtain the ratio of subparticle mass over progenitor par¬ 
ticle star formation rate and interpret this as a conversion 
timescale for each subparticle to become a star particle. We 
then randomly sample individual conversion times t for sub¬ 
particles using an exponential distribution of the appropri¬ 
ate timescale. If t < 100 Myr, a subparticle is deemed to be 
converted into a star particle with tage = ti, a mass equal 
to the subparticle mass and the Z* value inherited from the 
parent particle. 

For our analysis we use mdis = 10“^ Mq as our target re¬ 
sampling mass resolution. This represents a more reasonable 
Hll region mass, to better reproduce the discretisation in ob¬ 
served galaxies. To illustrate the effect of re-sampling, the 
left-hand panel of Figure [XT] shows the g — r colour and age 
of individual star particles in an actively star-forming EAGLE 
galaxy. We also plot the cumulative colour of the galaxy by 
star particle age as a black line. This shows the influence 
that the few youngest star particles can have on the colour 
of the entire galaxy. The left and right panels show the re¬ 
sult with and without re-sampling respectively, yielding dif¬ 
ferent overall colours for the galaxy (seen as the rightmost 
point of the black line). The presence of a single very young 
star particle in the simulation output causes this galaxy to 
appear ~ 0.2 mag bluer than with re-sampling. Though re¬ 
sampling can have a significant effect on individual galaxy 
colours, the colour distributions for the entire EAGLE pop¬ 
ulation are only marginally affected. At the low-mass end, 
the re-sampling generally serves to tighten the g — r colour 
distribution, move the blue peak to slightly bluer colours 
(~ 0.05 mag) and to suppress extremely blue outliers. 


APPENDIX B: COLOUR CONVERGENCE 

As simulation Recal-25 has a factor of 8 finer mass res¬ 
olution than the fiducial model Ref-100, the stellar mass 
threshold above which galaxies are considered well resolved 
is pushed to lower masses. By comparing colour distribu¬ 
tions of the Ref-100, Ref-25 and Recal-25 simulations (Ta¬ 
ble]^ S15) for galaxies within a certain mass range, we at¬ 
tempt to decouple the effects of simulation volume and reso¬ 
lution on the colours of low-mass galaxies in EAGLE. In Fig¬ 
ure |B^ we compare colour distributions for galaxies of mass 
9.45 < log(M./M 0 ) ^ 10.05 and 8.7 < log(M*/M 0 ) ^ 9.3 
in the top and bottom panels respectively. The histograms 
for differing simulation volumes have different y axis ranges, 
with the 25 Mpc simulation axis range a factor of 64 smaller 
to account for the differing simulation volumes. 
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Figure Al. Demonstration of the re-sampling technique for the stellar component of an example EAGLE galaxy from the 50^ Mpc® 
box at redshift z = 0.1. The left panel shows particle SSP g — r colours as a function of starburst age for an individual EAGLE galaxy 
as points coloured from blue to magenta indicating low to high metallicity. The black line indicates the cumulative galaxy colour when 
including particles with increasing age (from left to right) so that the far right point of the line gives the total galaxy colour. The right 
panel shows the same galaxy after star particles of age < 100 Myr have been re-sampled. These particles (black stars; only one in this 
example) are removed from the photometry calculation and the re-sampled stars (grey points) are included. This changes the cumulative 
g — r indicated by the black line, and yields a different overall colour. 


In the 9.45 < log(M*/M 0 ) ^ 10.05 mass range, the po¬ 
sition of the red and blue peaks appear roughly the same 
in the different simulations. However, the relative strengths 
of the red and blue populations differ, with the red se¬ 
quence being significantly weaker than the blue cloud in 
the high-resolution Recal-25 model compared to Ref-100 and 
Recal-25. This is consistent with the lower passive fractions 
in the high resolution simulation at a = 0.1 shown in S15. 

The 8.7 < log(Mt/M 0 ) ^ 9.3 range shows less con¬ 
sistency, with the red sequence becoming practically absent 
in the Recal-25 model while remaining in the Ref-100 and 
Ref-25 models. The redder colour and larger scatter of the 
blue population in the reference model is attributable to 
poor sampling of star forming gas in these galaxies. The 
lower star formation rates in the fiducial volume may also ac¬ 
count for the different colours. However, we also see a larger 
difference between Ref-25 and Ref-100 here, particularly in 
the relative contributions of the red and blue populations. 
We attribute the higher contribution of the red sequence in 
the Ref-100 model to the presence of large cluster environ¬ 
ments in the Ref-100 simulations, and thus quenched satel¬ 
lite galaxies, that are not sampled by the Ref-25 box. This 
suggests that volume effects also contribute to the weaker 
red sequence seen in the Recal-25 box. In both plots the 
greater area under the Recal-25 histogram is indicative of 
the systematic shift in galaxy number densities between the 
simulations, also seen in Figure 


APPENDIX C: SSP PARAMETER INFLUENCE 

As intrinsic galaxy colours are sensitive to star formation 
histories and elemental abundance patterns, comparing EA¬ 
GLE model colours directly to observed galaxy colours is a 
difficult way to disentangle the influence of different SSP 
parameters and to identify the source of any discrepancies. 


To go some way towards assessing how the EAGLE stellar 
metallicities and star formation histories influence our mock 
photometry, we use the simple photometric model without 
dust (N). Two sets of photometric data are first generated 
for the simulated galaxy sample using simulation output for 
one parameter while using empirical relations for the other. 
The galaxy metallicities and light-weighted ages (LWAs) as 


functions of stellar mass presented by Gallazzi et al. (2005 \ 


are used to provide the empirical input. For the LWA val¬ 
ues, we include a Gaussian scatter about the median values 
of the published width, which is assumed to be uncorrelated 
with metallicities. Glearly the assumption that galaxies may 
be treated as a single starburst and that the metallicity and 
age parameters are uncorrelated are poor, so the amount 
of information that can be drawn from this type analysis 
is limited. These plots serve as a basic qualitative illustra¬ 
tion of the influence of different SSP parameters on galaxy 
colours. 

The colour distributions in four M* bins are plotted for 
the simulation-empirical hybrid photometry models, and are 
compared to the eagle photometry in Figure [CT| The black 
lines indicate distributions of Ref-100 simulation galaxies. 
The distributions using the raw emission model with EAGLE 
ages and metallicities are plotted as solid histograms. The 
photometry models using observed LWA and Zt, values are 
plotted as the dashed and dotted lines, respectively. The 


observational data of Taylor et al. (20141 is also plotted in 
blue for comparison. 

We see that the age parameter has the biggest influence 
on the colour distribution, with the empirical ages introduc¬ 
ing a generally larger spread than metallicities, when com¬ 
pared to the pure EAGLE photometry. We have verified that 
this is still the case when we include the scatter on observed 
metallicity values. Figure |G1| shows that giving galaxies a 
single age stellar population using the observational LWA 


data of Gallazzi et al. (20051 (dashed line) works reasonably 
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data better than an empirical model assuming a single age 
population. 

The observational Z model reveals a poor fit to obser¬ 
vation for the two highest-mass bins. The red sequence is 
also much less prominent than seen in the observations and 
the other models across the M* range. 

The systematic effect of assuming uncorrelated scatter 
between the age parameter may also account for the fact 
that the colour distributions are broader and flatter than 
observed for this model, especially in the low-mass bins. 

The resolution effects that drive much improved agree¬ 
ment between observed low-mass colours and Recal-25 rel¬ 
ative to Ref-100 are noted in section]^ and appendix [B] In 
the lowest-mass bin of Figure [CT| we see that using observed 
metallicities has less impact on eagle colours than using ob¬ 
served LWAs. This indicates that star formation rate resolu¬ 
tion is the primary resolution effect on colours, with metal- 
licity resolution secondary to this. The presence of a faint red 
sequence is due to lower star formation rates and higher stel¬ 
lar ages than found in low-mass Ref-100 galaxies, whereas 
the position of the red sequence is redder by ~ 0.1 due to 
the higher than observed metallicities at these masses. The 
star formation rate resolution is also the main contributor 
to the redder than observed blue cloud position in Ref-100. 


Figure Bl. Comparison of g — r colour distributions of the 
Ref-100, Ref-25 and Recal-25 simulations (see table to as¬ 
sess the weak convergence and volume effects on model colours. 
Ref-100 and Ref-25 have the same resolution, while Recal-25 
has a resolution 8 times finer. The top and bottom panels show 
galaxy stellar mass ranges of 9.45 < log(M*/M0) Si 10.05 and 
8.7 < log(Mt/MQ) si 9.3, respectively. In both panels, the blue 
and red histograms represent the counts per colour bin in the 
100 Mpc and 25 Mpc simulations respectively. The Recal-25 and 
Ref-25 simulations, are plotted as solid and dashed lines, respec¬ 
tively. Separate y-axes are labelled and coloured to correspond to 
the two volumes, with their ranges scaled by a factor of 64 to ac¬ 
count for the differing box volume. Both resolution and box size 
appear to significantly effect the colour distributions of low-mass 
EAGLE galaxies. 


well in the two most massive bins where stellar populations 
are old. However in the lower-mass bins where galaxies are 
generally younger they provide a poor lit to the observed 
colours, inferior to our model photometry using the com¬ 
plex star formation histories of EAGLE (solid line). 

The bimodality seen for the full EAGLE photometry in 
the two most massive bins, but not for in the observational 
LWA model, shows that the EAGLE populations are intrinsi¬ 
cally bimodal in age. This supports the assertion that there 
is an excess of star forming galaxies in this regime relative 
to the observed population. The bluer than observed high- 
M* red sequence in the observational LWA model could be 
a result of the lower metallicities of high-M* galaxies. The 
inferior agreement of the observational LWA model relative 
to the full EAGLE model in the lower-mass bins suggests that 
the complex star formation histories of EAGLE reproduce the 
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Figure Cl. g — r colour distributions for EAGLE galaxies for the 
dust-free model. The 4 panels show colour distributions for 4 bins 
of stellar mass, as indicated by the legend. The solid, dotted and 
dashed lines show the EAGLE SSP values and the EAGLE SSP val¬ 
ues with Gallazzi et al. (2005) metallicities and ages respectively. 
Gallazzi metallicities and ages are assigned to each galaxy, based 
on the median parameter values at the galactic stellar mass. Z 
values are simply taken as the observed median value. The LWA 
values are sampled from a normal distribution with standard devi¬ 
ation taken from Gallazzi et. al. (2005), assuming that the scatter 
in age and metallicity is uncorrelated. We see that the complex 
star formation histories of EAGLE provide a better match to the 
observed colour distributions than a single SSP model using em¬ 
pirical values for age and metallicity. 
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